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We show that the stability of three-dimensional Skyrmions in trapped Bose-Einstein condensates 
depends critically on scattering lengths, atom numbers, trap rotation and trap anisotropy. In 
particular, for the 87 Rb |F = l,m/ = —1), |F = 2,m/ = 1) hyperfine states stability is sensitive to 
the scattering lengths at the 2% level, where the differences between them are crucial. In a cigar 
shaped trap, we find stable Skyrmions with slightly more than 2 x 10 6 atoms, a number which scales 
with the inverse square root of the trap frequency. These can be stabilized against drift out of the 
trap by laser pinning. 

PACS numbers: 03.75.Lm, 03.75.Mn 



I. INTRODUCTION 

A broad range of physics has been uncovered in the 
study of vortices in Bose-Einstein condensates (BECs) 
jH Q Q Therefore adding internal degrees of freedom [j 
to produce multi-component topological defects, such as 
coreless vortices and Skyrmions, should reveal further 
interesting phenomena. Skyrmions, which are the subject 
of this paper, were first introduced in particle physics 0, 
00 El where their application 10] has recently triggered 
a controversial search for penta-quarks [ll[ . 

It has been shown that Skyrmions can be energetically 
stable in a harmonically trapped two component BEC 
[T^ . and that the Skyrmion top ology c an be seeded by 
phase imprinting techniques [ijj, LLJ, LL3 • However their 
experimental realization remains to be achieved. In order 
to facilitate this, we have surveyed the relevant parame- 
ter space and identified regimes in which Skyrmions are 
stable. 

Other work on Skyrmions in BECs includes their study 
in a homogeneous background [l^ and in spin-1 spinor 
BECs . Trapped Skyrmions with winding numbers 
greater than one have also been shown to be stable [T^ . 

Skyrmions are topological solitons. Their two compo- 
nent wavefunction has a three dimensional structure that 
cannot be continuously deformed into a simpler struc- 
ture, such as a coreless vortex. A defining property of 
the Skyrmion is that the atomic field approaches a con- 
stant value at spatial infinity: for a trapped BEC this 
means sufficiently far from the trap centre. Then we can 
identify all points at infinity, turning the domain of our 
wave function into the compact space S 3 , the three di- 
mensional spherical surface in four dimensions. The two 
component wavefunction of the BEC at a particular posi- 
tion can be represented by the SU (2) matrix which maps 
a standard state, such as (1,0), into it. The wavefunc- 
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FIG. 1: (color online) Three dimensional density profiles for 
a Skyrmion with Ntot = 2 x 10 6 , N r = 0.6 in the presence of a 
pinning potential (Vb = 6.5huj x , wo = 1.3o~ x , see Eq. (J) and 
in a rotating trap (Q = 0.060;^). The anisotropy of the trap is 
cu±/cuz = 1.3 and the scattering lengths are the a%j. See text 
for definitions of symbols and for values of other parameters. 
The central (blue) torus is an isosurface of the line compo- 
nent. An isosurface of the ring component (red) is shown for 
x < 0: on the y-z plane between the isosurface sections the 
ring component density is indicated by a colormap from red 
(lowest) to purple (highest). The circulation associated with 
each vortex is also indicated. 



tion is then represented by a map from S 3 onto SU{2). 
A fundamental result of homotopy theory is that all such 
mappings can be classified by an integer topological in- 
variant, the winding number W []Jj- 

For a W = 1 Skyrmion, one component of the BEC 
carries a charge one line vortex, while the other com- 
ponent circulates around the first and through the line 
vortex core, carrying a charge one ring vortex singularity 
(Fig. EJ. The outer component carrying a ring singular- 
ity will be called the "ring component" (index 1), and 
the inner component with a line vortex will be called the 
"line component" (index 2). The filling of the ring vor- 
tex core by the line component stabilizes the dual vortex 
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structure as outlined in [T^ . 

Considering the specific case of 87 Rb hyperfine states 
|1) = \F = l,m/ = -1), |2) = \F = 2,m/ = 1), and 
identical spherical traps with frequency uo — 2tt x 7.8 Hz, 
we found a lower limit on the atom number for stability 
against cylindrically symmetric perturbations of about 
3 x 10 6 atoms in a spherical trap. This limit is slightly 
lower if the trap is elongated. Stabilization against gen- 
eral perturbations, which requires trap rotation, is only 
possible with at least 9 x 10 6 atoms. However with laser 
pinning, and slight adjustment of the scattering lengths, 
Skyrmions with as few as 2 x 10 6 atoms are found to 
be stable. The scaling relation given in the next section 
allows these results to be applied to traps with other fre- 
quencies. 

A result with particular significance for the experimen- 
tal realization of Skyrmions is that their stability depends 
on the atom-atom interaction strengths at the 2% level, 
where the differences between them are crucial. 

This paper is organized as follows: in section ITU we in- 
troduce the underlying physics and relevant parameters. 
Section IIII Al presents the results of our search for low 
atom number Skyrmions, section IIII Bl contains a brief 
discussion of trap offset between the two components, 
section ITvl shows how to stabilize Skyrmions against drift, 
and section gives an overview of their sensitivity to 
scattering lengths. The appendix discusses the origin of 
the trap offset between the components. 



II. METHOD 

To determine Skyrmion stability we numerically solve 
the Gross-Pitaevskii equation (GPE) in imaginary time, 
which corresponds to a numerical minimisation of the 
energy functional (see e.g. |2(J). The coupled GPEs for 
a two component BEC in a frame rotating with angular 
velocity ft compared to the lab frame are 



ftJ^iW = ( - |V + v(x) + u » 



3=1,2 



n-Lfi(x), 



(i) 



where for i = 1,2 denote the BEC wavefunctions 
of the two components, normalised to the number of 
atoms in each component, Ni and iV^. For later ref- 
erence we define N tot = Ni + N2 and N r = N2/N tot . 
V(jc) = m(cc;j_(x 2 + y 2 ) + cu 2 z 2 )/2 is the harmonic trap- 
ping potential, and m is the atomic mass. The interaction 
parameters Uij are proportional to the scattering lengths 
for intra- (U\\,U 22) and inter-component (E/12) interac- 
tions: Uij = 4:7rh 2 aij jm. L is the angular momentum 
operator. 

The imaginary time GPEs follow from Eq. by the 
replacement r = it. To solve them, we start with a trial 
wavefunction that is in the W = 1 topological class 12]. 
During the imaginary time evolution the normalisation 



of the wavefunctions is held fixed at Ni, N2. If the wave- 
functions converge, and are stable against general pertur- 
bations, the result is a local minimum of the energy. This 
is a steepest descent method in the function space with 
fixed norm, whereby the initial trial function is evolved 
in the direction where the energy decreases most rapidly. 

Experimental Skyrmion creation would involve phase 
imprinting a state in the W = 1 topological class onto 
the BEC jll El El- This would then relax towards the 
ground state of the class by means of thermal dissipa- 
tion. This process can be simulated numerically using a 
combined real and imaginary time formalism |2l|. 

After rescaling, x = x/a 2 , t = cj z t, ^ = ^iaj 2 j \[Wi, 
ft = il/u z , L = L/ft, where a z = yh/muo z is the har- 
monic oscillator scale length, Eq. has the dimension- 
less form: 



4*«<*)=( 



1 LO 



■-(x 2 + f) + z 2 



-n-L)*i(x), (2) 



.7 = 1,2 

where the interaction coefficients are defined as 



\ij = 4irNjdij \J moj z /h. 



(3) 



In a system with a spherical nonrotating trap the physics 
is fully determined by the Xij. All simulations used the 
parameters of 87 Rb. Early measurements gave the scat- 
tering lengths: an = 5.67 nm, CL22 = 5.34 nm, ayi = 
a,2\ = 5.5 nm |22| for the hyperfine states |1) and |2). 
After the bulk of our numerical work was done, we be- 
came aware of revised values for the scattering lengths: 
Ojm = 5.315 nm, a22 = 5.052 nm, a\2 — &21 = 5.191 nm 
|23l |24| . These differ from the a^- primarily by an over- 
all scale factor. Choosing this to be an /an = 1.067, 
according to Eq. (0 this is equivalent to an increase 
in atom numbers by the same factor, or by about 7%. 
The residual differences, although small, are significant. 
Therefore, our major results are presented for both sets 
of scattering lengths: early (a^ ) and revised (a^ ). 

In magnetic traps the two states will in general have 
slightly different trap frequencies due to the differential 
Zeeman effect HH23, which is reviewed in the Appendix. 
Taking gravity into account, this results in a small offset 
of the trap centres for each component which, for weak 
traps, can become comparable to the Skyrmion length 
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. We report in section IIII Bl that this 
effect should and can be minimized and therefore, un- 
less otherwise stated, we assume identical spherical traps 
with uj = uj z = cu± = 2tt x 7.8 Hz for both hyperfine com- 
ponents. In contrast, the change in non-linear interaction 
coefficients due to the actual different trap frequencies 
themselves is negligible. 

Although we present our results for specific trap fre- 
quencies, the scaling Eq. © implies that our results ap- 
ply for different trap frequencies provided the atom num- 
ber is scaled to keep the constant. For example if the 
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trap frequencies are increased by a factor of four, the 
atom number decreases by a factor of two. This scaling 
applies provided the Gross-Pitaevskii equation is valid, 
and hence fails if the density becomes too high. 

Skyrmions require that the interaction parameters Uij 
are such that the two components spatially separate. For 
a homogeneous system it can be shown that this corre- 
sponds to |2y, l2jj 



C/l 2 2 > U U U 2 2. 



(4) 



In our numerical survey of the BEC parameter space 
three dimensional spatial grids with 64 3 grid points were 
used. However the specific results presented in this paper 
were verified on bigger grids. To solve the GPE in imagi- 
nary time we used both the RK4IP algorithm [2^ and the 
corresponding first order interaction picture method. We 
found that the lower order method was faster. Typical 
imaginary timesteps were At = 10 -5 s. The simulations 
were done with the aid of the high level programming 
language XMDS [U . 



III. STABILITY AGAINST CYLINDRICALLY 
SYMMETRIC PERTURBATIONS 

A. Atom numbers 

Previous work had suggested that stable Skyrmions 
require total atom numbers of the order of 10 7 0. A 
major motivation for the work reported in this paper was 
to find stable Skyrmions with fewer atoms. We analyze 
stability in two stages. In the first stage, reported in 
this section, we consider stability against cylindrically 
symmetric perturbations. In particular, against collapse 
of the ring vortex, or deconfinement of the line compo- 
nent. In the latter case, the density of the line compo- 
nent drops off more slowly than the ring component at 
the edge of the condensate. Consequently the phase vari- 
ation associated with the line vortex extends to infinity 
and the wavefunction there cannot be identified with a 
single point in SU(2). Thus the topological classifica- 
tion of wavefunctions into classes with integer winding 
number fails. Numerically this is indicated by the calcu- 
lated winding number not being an integer. In the second 
stage we consider stability against non-cylindrically sym- 
metric perturbations. In particular, against drift of the 
line vortex singularity out of the trap. 

The results of our survey of stability against cylindri- 
cally symmetric perturbations are summarized in Fig. 
Points on the plots indicate whether a Skyrmion with 
the corresponding N r ,N tot values is stable (+) or not 
((g)). The two panels show the results for the two dif- 
ferent sets of scattering length introduced in section ITTI 
As expected from the discussion in that section, the re- 
vised scattering lengths require slightly more atoms for 
stability with given N r . The unstable cases with slightly 
too few atoms have imaginary time evolution for which 
the ring singularity (compare with Fig.|3J) moves inwards 
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FIG. 2: (color online) Skyrmion stability diagram against 
cylindrically symmetric perturbations, in a spherical trap: (+ 
blue) stable, ((g) red) unstable. The axes are the total atom 
number N to t = Ni + N2 and relative line component fraction 
N r = N2/Ntot- a) ciij set of scattering lengths. For N r > 0.6 
all tested cases were unstable, b) cnj set of scattering lengths. 
For N r > 0.65 all tested cases were unstable. 



until it reaches the Skyrmion core (defined as the volume 
with ring component flow through the line vortex core), 
and the flow velocity through the core increases. Even- 
tually, the flow velocity exceeds the local two component 
speed of sound (see Eq. ©) and the ring vortex collapses, 
destroying the ring cirulation. 

The filling of the vortex-ring by the line component 
can hinder the collapse of the ring vortex, and indeed 
we find stable Skyrmions with smaller Ntot if the relative 
fraction of the filling N r is increased, as shown in Fig. [21 
However, this works only up to N r ~ 0.65. For higher 
N r the line component escapes confinement by the ring 
component. Its phase variation means that the topolog- 
ical equivalence of R 3 and S 3 breaks down, and with it 
the integer valuedness of the winding number W. 

In summary, Fig. [3 shows a region with stability 
bounded towards lower N tot by ring vortex collapse, and 
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FIG. 3: (color online) Mach number of ring component in the 
y-z plane for a case near the stability boundary with Ntot — 
5 x 10 6 and N r = 0.5. The superimposed colormap shows the 
ring component density (blue - lowest density, red - highest 
density) . The peaks are at the location of the ring singularity. 
Scattering lengths ciij. 

towards higher N r by breakdown of the global topologi- 
cal structure. The lowest Skyrmion atom number allowed 
by these two constraints is N tot ~ 3 x 10 6 for the revised 
scattering lengths. The increase in central flow speed 
associated with the collapse of the ring vortex, in imagi- 
nary time evolution, led us to consider configurations in 
which this velocity was reduced. The flow speed is high 
through the core because it is narrow, its width being de- 
termined by the healing length, Fig. [31 The return flow 
speed of the ring component through the outer shell of 
the Skyrmion is typically orders of magnitude slower be- 
cause it occurs through a correspondingly greater cross 
sectional area. 

The quantisation of the circulation for a charge one 
ring vortex requires the velocity vi to satisfy 

/v r dl=-, (5) 

Jc m 

where C is any closed path around the ring singularity. 
Due to the high flow speeds through the core, most of 
the circulation accumulates there. Elongating the core, 
by stretching the trap into a cigar shape, should there- 
fore slow the core flow as the circulation is generated over 
a greater length of core. This is indeed what we found. 
However decreasing the atom numbers for which stable 
Skyrmions exist by stretching the trap reaches a limita- 
tion as the outer shell of the ring component thins and 
the line component pushes through to the outside. 

To prevent this we decreased CL22 by 2%. Achieving 
this might require scattering length manipulation, per- 
haps using an optical Feshbach resonance 30]. This in- 
creases the surface tension, as will be described in sec- 
tion and simultaneously makes the line component 
denser and more localized. With this we found a stable 
Skyrmion with N tot = 2.13 x 10 6 atoms and N r = 0.6 in 
a trap with anisotropy 77 = uj±/uj z = 1.3, uj±_ = 2tt x 7.8 
Hz. 

We also observed that increasing the fraction of line 
component 7V r , that has the stabilizing effect described 



above, lowers the central flow speed: Fig. 21 shows how 
the shape of the core varies for small deviations in N r 
and the effect of this variation on the flow speed. 

As long as the features of the Skyrmion are sufficiently 
larger than the wavelength of the sound-modes in ques- 
tion, we can approximate the background two component 
BEC as locally homogeneous and apply the standard re- 
sult for the speed of low energy excitations 4] 

c± = ( -— (Uiim + U 2 2n 2 



± VW"1 - ^2 2 n 2 ) 2 +4711712^) j . (6) 

Here ri\ and rt2 are densities of the ring and line com- 
ponents respectively. Demanding that both branches of 
the two-component speed of sound be real yields the 
condition JJyi < U\\U 22, which is violated in the ho- 
mogeneous Skyrmion system, Eq. (gj). We therefore 
only consider the real branch of Eq. ®. We define 
the Mach number Mi = Vi(x)/c+(x), as the flow ve- 
locity of component i expressed in units of the local 
two-component speed of sound. Fig. shows Mi along 
the Skyrmion core, the core width (defined as d(z) = 
{[J (x 2 +y 2 )\^i\ 2 dxdy]/[f l^i^dxdy}} 1 / 2 where the inte- 
gration is over the inner core area only), its density, and 
the unrescaled flow speed for a typical Skyrmion close to 
the stability boundary. The speed of the ring component 
though the thin vortex core of the line component has a 
sharp peak in the centre and reaches Mi ~ 0.15. 



B. Offset of the two trap centres 

As mentioned in section [H] and discussed in the Ap- 
pendix, the hyperfine states of the two-component BEC 
will in general experience differential gravitational sag, 
and hence slightly offset traps. We numerically investi- 
gated this by introducing a trap offset in the z-direction 
d = dz. We found for a Skyrmion far from the stability 
boundaries of Fig. b) (N tot = 9 x 10 6 and N r = 0.5 
[l2| with a radius of approximately 40 microns) that an 
offset of d = 2 jam produced instability. The hyperfine 
components spatially separate and the topological struc- 
ture breaks down. However, with an offset of d = 0.1/xm 
the Skyrmion remained stable. 

As the offset scales as uo~ 2 (see Appendix) it can always 
be made sufficiently small compared to the Skyrmion 

— 1/2 

length scale a z oc uo z by using a tighter trap. 
Skyrmions in tight traps also require lower atom num- 
bers, by the scaling argument given in section^ Another 
option is to operate at the magnetic bias field where the 
states experience the same trapping frequencies, which 
for our states is B = 3.23 G [23|. In the Appendix we re- 
late the trap offsets investigated above to magnetic field 
strengths in the vicinity of this value. Finally, it might be 
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FIG. 4: Influence of changes in iV r on the core (N to t = 5 x 
10 6 ): (Dashed) iV r = 0.47, (solid) N r = 0.5, (dot-dashed) 
N r = 0.58. Scattering lengths cnj. a) Shape of the core. 
Shown is the cross section of an isosurface of the density of 
the line component (712 = 2.5 x 10 -19 m -3 ) with the x-z plane, 
b) Mach number of the ring component Mi along the core for 
x=0, y=0. 
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FIG. 5: Skyrmion core properties (Ntot = 5 x 10 6 and N r — 
0.5). Scattering lengths cnj. a) Core width d defined in the 
text, b) Density of ring component, c) Flow velocity v\. d) 
Mach number Mi of the ring component along the core. 



advantageous to create Skyrmions in a tight optical trap, 
which would also widen the range of suitable hyperfine 
states to include \F = 2,m/ = 2) and \F = l,m/ = 1), 
between which exists an interspecies Feshbach resonance 

□ 



IV. STABILISATION AGAINST DRIFT 

Line vortices in single component trapped BECs are 
unstable towards drift out of the trap. However they may 
be stabilized by rotating the trap [l|, . The Skyrmion 
line vortex behaves in the same way [I3 • 

As there is no obvious experimental way to rotate the 
components independently, we rotate both. While the 
filled vortex in the line component is stabilized by rota- 
tion faster than a critical value Q c ,h there is a destabi- 
lizing effect on the toroidal flow in the ring: for rotation 
frequencies Q bigger than a critical value fi Cjr , the imag- 
inary time evolution creates an additional line vortex in 
the ring component. Due to the associated phase vari- 
ation the topological equivalence of R 3 and S 3 breaks 
down and with it the Skyrmion. 

Therefore stabilization requires Q Cj i < Q < ^ c , r . The 
critical frequencies Q Cj i and Q Cjr are determined numer- 
ically. To establish whether there is a stable Skyrmion 
for a given rotation frequency, we shift a previously ob- 
tained stable cylindrically symmetric state both paral- 
lel and perpendicular to the line vortex and evolve it in 
imaginary time. We observe either stabilization close to 
the original symmetry axis, or vortex drift out of the trap, 
or the disruption of the ring component circulation due to 
the creation of additional vortices. The response to trap 
rotation is studied for both sets of scattering lengths and 
varies significantly. The resulting windows of stable ro- 
tation frequencies are shown in Fig. with the shape of 
the data points distinguishing the different behaviour in 
imaginary time outlined above. The critical rotation fre- 
quency for stabilization of the line vortex increases with 
decreasing atom number, as is the case for an unfilled 
line vortex in a pancake trap , while the critical rota- 
tion frequency for disturbances in the ring component 
decreases slightly for smaller numbers. Therefore the 
window that allows the stabilization of Skyrmions by ro- 
tation closes around N tot = 8 x 10 6 for the early set of 
scattering lengths (a^) and around Ntot = 9 x 10 6 for 
the revised ones see Fig. 

For low atom numbers, Skyrmions can be stabilized by 
a laser pinning potential |33[ of the form |3J : 



V(x) 



w(z) 




exp [-2(x 2 +y 2 )/w(z) 2 



— , z % 



Wo, 



(7) 



where A is the wavelength of the laser and wo, Vo are 
constants that describe its width at the focal point and 
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FIG. 6: (color online) Windows of stable rotation frequency 
for various total atom numbers; N r =0.5 thoughout. 
(+ blue) stable, ((g) red) unstable against additional vortex 
creation in the ring component, ((} red) unstable against line 
vortex drift out of the trap. Main diagram: cnj set of scat- 
tering lengths. Inset: dij set of scattering lengths. The grey 
shaded area indicates the location of the stability region for 



intensity, respectively. Using the early scattering lengths, 
we found a stable Skyrmion with N tot = 2 x 10 6 and 
a pinning potential of Vb = 10Hlu x , wo = l.Sa x (a x = 
(h/mujx) 1 / 2 ) and A = 450 nm. 

However, the pinning laser beam can also have a desta- 
bilizing effect on the Skyrmions: it reduces the density of 
the ring component in the core, increasing the flow veloc- 
ity, which can lead to collapse of the ring vortex. If laser 
pinning and rotation are used together, it appears to be 
possible to limit the disrupting effects of both methods. 
Fig. [7| shows an example, the same as in Fig.QI with the 
pinning amplitude reduced to Vb = 6.5Hlu x , and trap ro- 
tation of = 0.06^. The scattering lengths here are also 
the early values, d^. We verified that this case is stable 
for independent perturbation of any one of the parame- 
ters mentioned above by 3%. The revised scattering 
lengths also allow a small stable Skyrmion with com- 
bined rotation and pinning, but require N to t = 2.13 x 10 6 
atoms. 



V. SCATTERING LENGTHS AND FESHBACH 
RESONANCES 

Due to potential uncertainties in the 8T Rb (revised) 
scattering lengths, we investigated the effect of changing 
them on Skyrmion stability. For a case with 9 million 
atoms we found that variations of an by 5% resulted in 
Skyrmion instability. For increased an this was due to 
disruption of phase separation. For decreased an the line 
component moved outward through the ring component. 
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FIG. 7: (color online) Small Skyrmion with N tot = 2x 10 6 and 
N r — 0.6, stabilized against drift with rotation and pinning. 
Scattering lengths dij . a) Colormap of line component density 
(blue - lowest density, red - highest density) in the x, z plane, 
b) Ring component density, c) Density profiles of the line 
(dashed) and ring components (solid) along the x-axis (y=0, 
z—0) . Parameters as for Fig. Q 



This is because the relative magnitude of an and a 2 2 de- 
termines which component surrounds the other. In the 
Skyrmion configuration an is larger and hence lower den- 
sities of the ring component are energetically favorable. 
We did not find that increases in an were destabilizing. 
However, decreases of 2% caused instability, again due to 
loss of phase-separation. Similarly, decreases in a22 were 
benign, but increases of 2% lead to instability. To verify 
some of these results we used numerical grids with 256 3 
points. 

The sensitivity of Skyrmion stability to the scattering 
lengths can be understood in terms of the phase separa- 
tion criterion or more generally in terms of the surface 
tension a at the component interface |23, 133- For an 
interface within a homogeneous background, and in the 
limit of Uu = U22, |35| a is 
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U 12 



2m \ y/UuUm 
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The derivative of this with respect to Uu has a pole 
where a has a zero (the critical value for phase separa- 
tion). It then predicts large variations in surface tension 
for small variations in Uu. Even though we find signifi- 
cant deviations from the homogeneous results, like stable 
Skyrmions in the regime U 2 2 < U\\U 22^ this qualitative 
prediction from the homogeneous case holds. The zero of 
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FIG. 8: (color online) Skyrmion with Ntot 
N r = 0.5 for different surface tensions a. 
(solid black) normal a: a~ij set of scattering lengths, 
(dashed red) low a: an increased by 2%. 
(dot-dashed blue) high a: a±2 increased by 5%. 
a) Densities of the line- and ring components across the core 
(z=0, y=0). b) Mach number (Mi = vi/c+) of the ring com- 
ponent along the core (x=0, y=0). 



(7 is shifted from C/f 2 = U11U22 in the trapped case, but 
our simulations show big changes in surface tension for 
small changes in the Uij, as displayed in Fig. [5] Fig. 
a) illustrates the variation in the size of overlap regions 
between the line and ring components. Fig. b) shows 
the Mach number. The qualitative differences result from 
changes in the shape of the core. 



VI. CONCLUSION 

Some useful conclusions may be drawn from our study 
regarding the experimental prospects for producing sta- 
ble Skyrmions in trapped BECs. The situation is more 
optimistic than initial work might have suggested [T^ . 
as we have shown that substantially less atoms are re- 
quired for stable Skyrmions to exist. However, we have 
also found that the windows of stability can be rather 
small, and that they are sensitive to the precise values 



of the scattering length ratios. The ability to manip- 
ulate the scattering lengths, for example using optical 
Feshbach resonances, would improve the robustness and 
flexibility of experiments. Stabilization techniques such 
as laser pinning are also likely to be useful. 
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APPENDIX A: APPLICATION OF THE 
BREIT-RABI EQUATION 

This appendix outlines how the magnetic field depen- 
dence of the magnetic moments of the hyperfine states 
constituting the Skyrmion leads to different trap fre- 
quencies. The energy of an atom in the hyperfine state 
\F,rriF)i with electronic spin 1/2, in an external magnetic 
field B is given by the Breit-Rabi equation |25|, l3q 



E(F,m F ,B) 



± 



2(27+1) 



gHi B m F B 



1 



2m p 2 

^ QC I QC 



1/2 



B. 



(Al) 



The nuclear spin for 87 Rb is I = 3/2. Here A^f is the 
hyperfine splitting at zero magnetic field, gi and gj are 
the gyro magnetic ratios of electron and nucleus respec- 
tively and hb is the Bohr magneton, x is a dimension- 
less magnetic field and the upper sign is for states with 
F = I + 1 /2. the lower one for states with F — I — 1/2. 
For 87 Rb 0,111 

A hf = 6834.68261090434(3)MHz, (A2) 

gi = 0.9951414(10) x 10~ 3 , 

gj = 2.00233113(20). (A3) 

The following discussion is one dimensional for simplicity. 
We assume the magnetic trap has the form 



B = Bn 



B 2 x 2 



(A4) 
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where Bq is the bias field. We can expand Eq. (|A1|) 
around some Bq as 

E(i, B) = E^(Bo) + E?\Bo)(B - B ) (A5) 

+ E^\B )(B-B ) 2 + O(B-B ) 3 , 

where all information concerning the state has been gath- 
ered in the collective index i. Inserting the magnetic field, 
we obtain the following potential for atoms in state i 



V(x)W = Eg* (B ) + E^(B Q )B 2 x z + 0(x 4 ). (A6) 

The ratio of trap frequencies for atoms in states i and j 
at bias field Bq is therefore 



k 

2 

(j) 



E[ j) (B ) 



(A7) 



A gravitational acceleration g will displace the equilib- 
rium position in the trap from x = by = g/uJ^y 
The difference in the gravitational sag between the two 
hyperfine states is then 



\x, 



U) 



x, 



(i)l 



(j) 



2 



10 



(i) 



E[ j) (B ) 



(A8) 



If we apply this result to the 7.8 Hz trap and hyper- 
fine states considered in this paper, the resulting offset is 
about 11/xm at a bias field of 1 Gauss. Eq. (|A1|) predicts 
zero offset at B — 3.22946 Gauss. The offset is about 
2/mi at B = 2.85 Gauss and about 0.1/im at B = 3.21 
Gauss, which are the values considered in section ITlI Bl 
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